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Abstract — This paper presents a method for imaging of moving 
targets using multi-static SAR by treating the problem as one 
of spatial reflectivity signal inversion over an overcomplete 
dictionary of target velocities. Since SAR sensor returns can 
be related to the spatial frequency domain projections of the 
scattering field, we exploit insights from compressed sensing 
theory to show that moving targets can be effectively imaged 
with transmitters and receivers randomly dispersed in a multi- 
static geometry within a narrow forward cone around the scene 
of interest. Existing approaches to dealing with moving targets 
in SAR solve a coupled non-linear problem of target scattering 
and motion estimation typically through matched filtering. In 
contrast, by using an overcomplete dictionary approach we 
effectively linearize the forward model and solve the moving 
target problem as a larger, unified regularized inversion problem 
subject to sparsity constraints. 

Index Terms — Multi-static SAR, imaging, regularization, spar- 
sity 

I. Introduction 

Synthetic aperture radar (SAR) is a remote sensing system 
capable of producing high-resolution imagery of target scenes 
independent of time of day, distance, and weather Conven- 
tional SAR radars are monostatic, with collocated transmit 
and receive antenna elements. These SAR sensors coherently 
process multiple, sequential observations of a scene under 
the assumption the scene is static. When the scene changes 
between these observations, as occurs when objects move, and 
these changes are ignored, blurring, defocus, and other artifacts 
are introduced into the reconstructed imagery. This is because 
the Doppler shift of moving objects are then determined not 
only by their geometric location but also by their velocity. 
Imaging of scenes with moving targets has gained increasing 
interest as the desire for persistent and urban sensing has 
grown. 

Moving target localization has proven challenging in the 
case of single antenna conventional narrow-angle SAR uti- 
lizing conventional reconstruction methods, such as the polar 
format and the filtered-back-projection algorithms fH], [25, due 
to an inherent ambiguity in target geolocation and velocity. 
Consequentially, most techniques for imaging moving targets 
with conventional SAR aim at focusing and detecting smeared 
targets in SAR imagery lH), ID, O, Hi, 0. In recent 
years, however, a number of techniques to handle moving 
objects explicitly have been developed. Space-time adaptive 
processing (STAP) |6| exploits multiple-phase center antennas 
to suppress clutter and produce a moving target indication 



image. Velocity synthetic aperture radar (VSAR) |T| exploits 
the velocity information contained in phases of a sequence 
of images formed at multiple receive antenna elements. Dual- 
speed SAR |8| has the radar platform move sequentially at 
two different velocities during the radar data collection time. 
Distributed antenna radars also have the potential to break 
the velocity-location ambiguity due to multiple phase centers 
of the antenna, while at the same time substituting spatial 
diversity for conventional bandwidth resources. Recent work 
on multi-static |9| and the related MIMO [T0"|, fTTl, lEl 
radar with coherent processing has shown the potential for 
resolution improvement that can far exceed the limit suggested 
by conventional arguments based on the radar's waveform. 

The conventional approach to resolving the moving target 
localization problem is to perform matched filter reconstruc- 
tion at every pixel for every possible velocity hypothesis 
independently, yielding a large space-velocity cube |9|, |4|. 
A target is then placed at the locations of maximal energy 
focus in the space-velocity cube. An approach based on the 
inversion of the forward operator is presented in ifTSll . where a 
filltered-backprojection approach to imaging of moving targets 
takes the form of a weighted matched filter These approaches 
require the solution of many large inversions and result in a 
large, somewhat ambiguous output. 

In this paper we utilize an inverse problem formulation 
and insights from sparse signal representation and compressed 
sensing for effective imaging of dynamic environments us- 
ing distributed antenna SAR sensor geometries. The non- 
linear problem of the coupled target localization and velocity 
estimation is linearized by construction of an overcomplete 
dictionary of velocity states. The resulting inversion leads 
to a non-convex optimization problem which we efficiently 
solve through convex relaxation. In contrast to the filtered- 
backprojection and the matched filtering approaches of fl^ 
and [9|, in our approach, velocity estimation is performed 
during the image formation process and all velocity hypothesis 
are evaluated jointly in a single optimization framework. 

II. Observation model 

We consider a multi-static system consisting of widely 
separated transmit and receive elements within a forward cone 
positioned at the center of a scene of interest. We assume 
that different transmitters send out probes in TDMA fashion, 
while the receive unit coherently processes signals of all 
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receivers across all snapshots, i.e. pulses. The scene of interest 
is modeled by a set of point scatterers reflecting impinging 
electromagnetic waves isotropically to all receivers within the 
forward cone, thus, allowing for the coherent processing of 
all received signals. The reflection coefficient of the point 
scatterer is a complex number with an unknown amplitude 
and a random phase lfT4l . 

We introduce a coordinate system with the origin in the 
center of the area of interest and, for simplicity, model the 
scene as two dimensional. Fig. [T] illustrates this set up. The 
relative size of the scene is assumed to be small compared 
to distances from the origin of the coordinate system to all 
transmitter and receivers, such that transmit and receive angles 
would change negligibly if the coordinate origin moved to any 
point in the scene. Furthermore, we neglect signal propagation 
attenuation. 

^ Y 

^-th transmkfer 




Fig. 1. Geometry of the fci-th transmit-receive pair witli respect to tlie scene 
of interest. All transmit and receive pairs are restricted to lie within a forward 
cone of the angular extent AS. 



A. Stationary scene model 

The complex signal received by the /-th receiver for the 
excitation from the fc-th transmitter reflected from a point 
scatterer at the spatial location x = \x^yY is given by 

"TkliPl = s(x) 7fe {t - rfc;(Xo) - Tfc;(x)) , 

where s(x) is the reflectivity of the scatterer, jk{t) is the trans- 
mitted waveform from the fc-th transmitter, Xq ~ [0, 0]^ is the 
scene's origin and Tki{xo) is the sum of the signal propagation 
delay from the fc-th transmitter to the scene's origin, Tfc(xo), 
and the propagation delay from the scene's origin to the l-th 
receiver, t;(xo), so that rfe,(xo) = Tfe(xo) +Ti(xo). Under the 
far field assumption, the propagation delay from the scene's 
origin to the scatterer at the location x is given by: 

rfci(x) = Tfc(x) + T((x) = --x^{ek + e;) = --x^eu, (1) 

c c 

where e^j = [cos^fc, sin0j.]^ and e/ — [cos 0/, sin <^;]-^ are 
unit vectors in the direction of the fc-th transmitter and Z-th 
receiver respectively. The signal delay to the location x in 
the direction of the fc-th transmitter is rfe(x) = —^x^ek and 



similarly, the signal delay to the location x in the direction of 
the l-th receiver is r;(x) = — ix^e;. Thus, the propagation 
delay is determined by the projection of a scatterer's location 
onto the kl-th transmit-receiver pair's bi-static range vector 

For extended scenes, multiple scatterers will have the same 
projection onto the bistatic range direction. The collection of 
such scatterers satisfies {xjix-'^efc; = p} . These scatterers are 
simultaneously illuminated and have their collective response 
Qki (p) registered at the receive antenna with the same delay. 
The so-called range profile qki (p) is an aggregate response at 
each delay or range and is given by: 

Qklip) = I s(x)(5 ( p - ix'^Gfc; )dx. 

Jllx|l<L V ^ / 

Under the far field and narrow-band transmit signal assump- 
tion ifTSll . the overall received signal from the entire ground 
patch is assumed to be a superposition of the returns from all 
the scattering centers and is given by: 

rki{t) = J Qkiiphk (^t - Tki{Xo) + dp. (2) 

In terms of the spatial reflectivity function s(x) the received 
signal is given by: 

rki{t) = / s(x)7fe (t - rfci(xo) - rfci(x)) dx, (3) 

J ||x||<L 

where Tf^i{x.) is given in ([T]). 

B. Moving scene model 

When moving scatterers are present in the scene we are 
interested in producing a focused image of the spatial re- 
flectivity function at some reference time Uef- We assume 
that the scatterer at location x has an associated arbitrary 
constant velocity vector Vx = [vx,Vy]'^. Note that in the 
case of non-constant motion, the true scatter motion can be 
well approximated as constant when the time scale of the 
coherent processing interval (CPI) is relatively small. We 
consider the general case when the CPI interval contains 
multiple transmitted pulses. 

Let us first consider the effect of the scatterer motion during 
one pulse transmission following the analysis developed in 
1 1 3 1 . In the case of motion, the delay of the transmitted 
waveform is dependent on both the location and the velocity 
of the scatterer The signal reflected from a single moving 
scatterer will now be of the form 

rkl{t) = s(x)7A; (t - Tfc;(Xo,X, Vx)) , 

where Tfe((xo,x, Vx) is the delay including the effects of 
motion. Next, we derive an expression for t — Tfc;(xo,x, Vx) 
in the case of motion. 

Let ix.vx(0 denote the time when the transmitted wave 
reaching the l-th receiver at time t interacted with the scatterer, 
that at the reference time Tfc(Xo) was located at x. Recall that 
Tfc(xo) is the propagation delay between the fc-th transmitter 
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and the scene's origin. Thus, the scatterer at time tx,vx(i) is 
actually located at x + Vxtx,vx(i)- We can write 



Tz(Xo) 



MY 



(4) 



Solving the above equation for tx,vx(i) and noting 
that Tfc(x;tx^v^(0) = --J [x + Vxtx,v,(i)] and 
T;(x;?x,v^(0) = -Hx + Vxix,v^(t)] e/, the ai-gument 
of the delayed pulse is derived to be; 

t - Tfe;(Xo,X, Vx) = 

= i - [Tfe(xo) + Tfc(x;fx,v^(i))] 

- [Ti(Xo) +T;(x;ix,v,(i))] 
= -rA;(xo) + x'^efc/c 

+ - {t - Tii^o) + x^e,/c 

1 - v^e;/c 

w t - Tki{yio) - m{y^) 

T / T 

, Vx^ efc( / x^ 
H i - Tfc(Xo) H 



where we used the fact that |vx|/c ^ 1, such that 
, ^ ; « 1 + Vx efc;/c. 

Extending the model to the case of multiple probe transmis- 
sions during the coherent processing interval (CPI), we assume 
that the scene is imaged at some arbitrary time tref outside 
the given pulse interval. Then, the scatterer located at x at the 
reference time tref, will be located at x + Vx(<a: — tref) at 
time tk that represents the time when the /c-th transmitted pulse 
reached the scene origin. Updating x of the previous equation 
with X + (ife — ire/)vx, we obtain for the fc-th transmitter 
pulse: 



t - Tfe;(Xo,X, Vx) « 

t - Tki{yio) - Tkl{-X.) 



(5) 



T 

Vx eki 



{t + tk- tref + efci(x,Vx)) 



where efci(x,Vx) = -Tfe(xo) + [x + {tk - ire/)vx]^ei/c. 

Finally, the forward observation model in the presence of 
motion becomes: 



rki{t) 



|x||<L 



S(x)7fe {t - Tfe;(Xo,X,Vx))dx, 



(6) 



with the argument t — rfc;(xo,x, Vx) given in ^ and s(x) 
representing the spatial reflectivity function at the reference 
time. Comparing this equation to the received signal model 
for the stationary scene given by (O, we see that the two 
models differ by an additional delay attributed to the scatterers' 
motion. The additional delay of a scatterer present at the 
location x at the reference time tref is proportional to the 
projection of the scatterers velocity Vx to the fc^-th transmit- 
receive pair bistatic range vector eki and the time interval in 
between the observation time and the reference time. 

For narrowband waveforms, defined by 7fc(t) — 
lk{t)e~^'^''* , where 7fc(i) is the low-pass equivalent signal 



and Ldk the carrier frequency, (|6]l is approximated |T3l by: 

rki{t) 



s(x)7fc (t - Tfc;(Xo) - Tfc;(x)) 



||x||<L 



-Tfci(x) + [t+tfc-*re/+efcl(x,Vj;)]- 



dx. (7) 

The additional phase shift of (|7| is a function of the quantity 

t^k- 



'-, which is basically the Doppler shift. Thus, the 
received signal is in the familiar form of the superposition 
of time-delayed and Doppler-shifted replicas. Notice that the 
Doppler shift is unique for each transmit-receive pair as it 
depends on the bistatic range vector eki- 

C. Discrete model 

A discrete version of the model in (|6]l or (|7]i can be obtained 
by discretizing the spatial variable x and sampling in time 
which, in the presence of receiver noise n, becomes: 

p 

r = ^*p(Vp)sp + n = *(V)s + n. (8) 

p=i 

In this equation, r represents the observed, thus known, set 
of return signals at all receivers across time. Its elements are 
indexed by the tuple (fc, I, tg), with tg being the sampling times 
associated with the fc^-th transmit-receive pair. The reflectivity 
of the p-th spatial cell or pixel is denoted by Sp G C^^^ and 
#p(vp) is the vector capturing the contribution to the received 
signal of a reflector that was located in the p-th pixel at the 
reference time tref and moved with the constant velocity Vp 
throughout the coherent processing interval. Stationary point 
reflectors are included in this model by simply setting Vp = 0. 

The received signal model described in (jSj, represents the 
observation model of unknown scatterers' reflectivity coeffi- 
cients Sp and their corresponding velocity vectors Vp. While 
the scattering coefficients Sp enter the problem linearly, the 
unknown velocities Vp do not, so the overall problem is 
nonlinear and coupled. When the velocities are known the 
remaining equation for Sp is linear, however, and straight 
forward focusing and estimation of scattering coefficients is 
possible. When the velocity is ignored (set to zero) or set 
to an incorrect value, the resulting reconstruction exhibits 
defocusing of the energy of the moving scatterer |[T], I?). These 
observations lead to one conventional approach to solving this 
problem. In particular, reconstructions are performed for every 
possible velocity yielding a large space-velocity cube 19], ||2l, 
[4]. Velocity slices where the image is well focused are as- 
sumed to indicate the correct velocity at a pixel. This approach 
requires the decoupled solution of many inversion problems 
and results in a large, somewhat ambiguous output. Recently, 
an approach based solely on the inversion of the forward 
operator is presented in liT3l . where a filtered-backprojection 
approach to imaging of moving targets takes the form of a 
weighted matched filtering which still produces a large space- 
velocity cube. In this paper we present a different approach 
based on recent results in sparse signal representation using 
overcomplete dictionaries ifTSl , sparsity based reconstruction, 
and compressed sensing lfT6l . ifTTl . which is described next. 
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III. OVERCOMPLETE DICTIONARY APPROACH 

Sparse signal representation aims at capturing a complicated 
signal as a linear combination of a few generating elements 
|IT5|. In particular, a signal in a given class of dimension M 
should be representable by a small subset of a collection of 
Q generating elements. For the case when M = Q and the 
elements are independent, the collection is termed a basis. 
When Q > M the collection is termed an overcomplete basis 
or dictionary. To represent this problem mathematically, let 
r e TZ'^'^^ represent the signal, 4> e TZ^^^ the dictionary 
and s £ TZ'^^^ represent the linear coefficients, such that r = 
$s. Since # has a non-zero null space many solutions are 
possible. What is sought is a sparse solution with only a few 
non-zero elements. 

While optimal design of dictionaries $ is a topic of general 
interest, in this work we assume the dictionary is given and 
fixed based on prior knowledge of the expected velocities in a 
scene. The problem is then one of finding an optimally sparse 
solution. A direct formulation of this problem can be given 
as: 

min ||s||o s.t. r = $s, 

s 

where || • ||o denotes the Iq norm, which counts the number 
of non-zero elements of the argument. Unfortunately, this 
formulation is computationally difficult to solve, as it involves 
NP-hard enumerative search which is prohibitively expensive 
for even moderate sizes of Q. A number of alternative, indirect 
techniques have been developed, based either on relaxation 
techniques or iterative greedy algorithms. 

The convex relaxation approach relies on the fact that 
besides the norm, the li norm also promotes sparsity in 
a solution. The li norm is defined as ||s||i = ^X^E^i ' 
where (s)i is the i-th element of s. This norm is a convex 
function of its arguments. The relaxed version of the problem 
then takes the form: 

min ||s||i s.t. r = $s, 

s 

which is essentially a linear program (LP). The use of this 
formulation has also been motivated by the fact that under 
certain conditions on the overcomplete basis the original 
problem and the relaxed version can be shown to have the 
same solution ifTSl . 

When the signal r is noisy, the signal representation problem 
becomes a signal approximation problem. The convex relax- 
ation formulation of the noisy signal approximation problem 
is given by: 

min ||s||i s.t ||r - *s||2 < 5 

s 

where 5 represents a small noise allowance. Instead of satis- 
fying the relationship exactly, the solution coefficient vector 
s is allowed to satisfy the relationship approximately. This 
problem is known in the literature as noisy basis pursuit |fT9l . 
We make use of this formulation in what follows. 

Compressed sensing (CS) ||T6| . ifTTI . takes the sparse rep- 
resentation framework one step further by seeking to acquire 
as few measurements as possible about a sparse signal, and 



given these measurements, reconstruct the sparse signal either 
exactly or with provably small probability of error using 
formulations like the above. Most of the work in CS assumes 
that the projections are drawn at random. However it is also 
known that Fourier measurements represent good projections 
for compressed sensing of sparse point like signals [l?]. This 
result immediately connects to radar measurements, as SAR 
sensors can be viewed as measuring samples of the stationary 
scattering field in the spatial frequency domain. 

A. The new formulation 

We exploit the overcomplete dictionary approach to signal 
representation to create a new formulation of the dynamic 
SAR inversion problem. In particular, we first introduce an 
appropriate overcomplete dictionary of velocity hypotheses 
by constructing a grid of all possible scatterer velocities at 
each location. This defines an over-complete representation of 
the received signal r. This representation is combined with 
the SAR observation equation resulting in an inverse problem 
that is linear in terms of an extended reflectivity coefficient 
vector, but is now under-determined. We then apply a modified 
version of the sparsity seeking formulation of (fTTI) to solve the 
resulting large, under-determined linear problem. Details are 
given next. 

Step 1: Dictionary Definition. We hypothesize that a scat- 
terer velocity VpjVp belongs to one of a discrete set of 
velocities V: 

Vp e V = {Vl = 0,V2, V3, ■ ■ • , Vat}, 

where N denotes the size of the velocity grid. The first 
velocity vector is set to zero to allow for stationary targets. 
The original observation vector *p(vp) describing the con- 
tribution of the p-th pixel to the received signal r in terms 
of an unknown pixel velocity Vp now becomes the matrix of 
[$p(vi), $p(v2), . • • *p(v7v)] composed of the contribution 
of each possible known velocity hypothesis at pixel p to the 
received signal. There are no unknowns in this matrix, in 
contrast to the original observation vector 

Finally, by combining this model at each pixel we obtain 
an overall overcomplete forward operator: 

*v = [*i(vi),...,*i(vjv),.-.,*p(vi),...,*p(vjv)]. 

To match this change, the reflectivity coefficient at pixel p 
is accordingly expanded to become a reflectivity coefficient 
vector: 

bpi 



ppN_ 

where the auxiliary variables 6p„ are constrained to the binary 
set hpn G {0, 1} to represent a 'true' or 'false' hypothesis that 
the p-th spatial location moves with n-th quantized velocity 
vector v„. Additionally, the variables 5p„ need to satisfy 
X^n ^pn ^ 1 for ^^ch p to recover the model of ([8]). This 
additional constraint specifies that at most one spatial reflector 
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is present within the p-th resolution cell at the reference time. 
Other, more complicated, models are also possible. 

Stacking up these single-pixel reflectivity coefficients yields 
an overall, extended reflectivity vector for the entire image: 



Sb = 



Combining this extended reflectivity vector with the overcom- 
plete forward operator yields our new overall linear motion- 
SAR forward model: 



r = *vSb 



n + q, 



(9) 



where Sb is the extended reflectivity coefficient and q rep- 
resents an additional noise term due to the velocity space 
quantization. In this new model is completely specified 
and contains no unknowns. We have essentially converted the 
original, difficult, non-linear problem to a selection problem 
with a linear observation. 

Step 2: Image Formation. We now treat the problem of 
image formation as a problem of inverting (|9]l. In principle, we 
seek to find a sparse vector Sb that best describes the received 
signal. In particular, we seek a solution of the following 
optimization problem: 



mm 

Sb ,b,s 



Sb 



subject to ||r - *^Sb||2 < <5 

(Sb)(p-l)Ar+n 

^6p„ = l p=l,...,P 

n 

6p„e{0, 1} n^l,...N. 

where 5 represents some noise allowance aimed at capturing 
both receiver and quantization noise and the Zo norm || • ||o 
counts the number of nonzero components in its argument. 
Since the coherence in typical SAR scenes is contained in the 
field magnitude (due to random phase), the Iq norm is applied 
explicitly to magnitudes of the extended complex reflectivity 
field |sb|. 

B. Solving the formulation 

The above optimization problem is a non-convex, mixed- 
integer program, known to be NP hard to solve. Thus, we 
resort to a two step procedure based on convex relaxation for 
its approximate solution. 

Step l:We first solve the following relaxed optimization 
problem: 



Sb 



arg mm 

Sb 



Sb 



subject to ||r — ^'^Sb||2 < e 



where |||sb|||i = Eti ^m^y^W + (^K).)' 
and (sb)i is the i-th element of the vector Sb. Note 
in this step we use the convex li penalty and do 
not enforce the exclusivity of velocity hypotheses 
represented by the auxiliary binary variables 6p„. 



Step 2:We then find the final reflectivity and velocity vec- 
tor estimates by selecting the maximum magnitude 
response at the p-th pixel: 



(sp,n) 



argmax|(Sb)(p_i)Ar+„| 



(10) 



This is consistent with the assumption that only one 
moving reflector is present in the p-th pixel at the 
reference time. 
By constructing the overcomplete dictionary and applying 
convex relaxation we have linearized the forward model at 
the expense of an increase in the size of the problem. The 
minimization problem in Step 1 is a second-order cone pro- 
gram, which we solve by a specialized large scale interior- 
point method for complex variables proposed first in f20l. 
The method is a specialized central path interior-point method 
with an approximate search direction found through a pre- 
conditioned conjugate gradient method which yields efficient 
solution of such large problems. 

Our approach relies on two levels of sparsity. First, the intro- 
duction and use of an overcomplete dictionary mandates that 
we seek a sparse solution vector Sb. In addition, compressed 
sensing theory relates the number of measurements necessary 
for accurate recovery of Sb to its underlying sparsity |16l, 
1 17], indicating that it may be possible to perform accurate 
recovery of sparse scenes with relatively few transmitters and 
receivers. 

We also want to emphasize that as long as the scene contains 
a sparse set of reflectors, there is no constraint on the object 
velocity the method can handle. The velocity grid resolution 
should simply be matched to the coherent processing interval 
and the maximal carrier frequency of the transmitted signal 
in order to avoid phase wrapping of the exponent in (|7]i. The 
phase error due to the velocity quantization is linear in Vx, 
so the velocity grid need not be constant: at smaller velocity 
magnitudes it can be coarser and at higher velocity magnitudes 
it can be finer Finally, the discrete model of dHJ implies that 
only a single scatterer is present at the reference time in each 
spatial pixel. For the case when the spatial grid is sufficiently 
coarse to include multiple point reflectors within the resolution 
cell, ^ can be written as: 



M[p) P 

E 

m — 1 p—1 



(11) 



where Vp „i models the velocity of the m-th reflector in the 
p-th resolution cell, that contains M{p) scatterers out of which 
one is allowed to be stationary. The algorithm accommodates 
the new model with a small change. In case of M(p) reflectors 
per spatial grid cell, the constraint &p„ = 1, should be 
replaced with X^n^P" ~ ^I{p)- After convex relaxation, 
([Tol l should be accordingly modified to retain M{p) maximal 
magnitude reflectivity values and their corresponding velocity 
vectors. 

IV. Transmitted waveforms 

We presented the overcomplete dictionary reconstruction 
algorithm without specifying the transmitted signals. The 
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formulation we presented is quite general and can potentially 
work for many waveforms and sensor configurations. 

In this initial work, to demonstrate the algorithm we con- 
sider multiple snapshot, non-overlaping transmissions of chirp 
and ultra-narrow band signals. Such signals are known to 
lead to a simple Fourier relationship between the stationary 
reflectivity field and the measurements for both monostatic and 
multiple distributed antenna configurations fTTi, f9l, (21). The 
corresponding forward operators posses compressed sensing 
properties, allowing for good reconstruction of sparse fields 
with few measurements (22], (TT\. In the following, we discuss 
the specific forward radar model resulting from transmission 
of such waveforms. 

The chirp signal is the most common spotlight SAR pulse 
lfr4J . given by 



' 2 — '' — 2 







Otherwise, 



where ojk is the center frequency and is the so-called chirp 
rate of the k-th transmit element. The frequencies encoded by 
the chirp signal extend from ujk — o-kTc to LOk + akTc, such that 
the bandwidth of this signal is given by Bk — Ultra- 
narrow band waveforms are special cases of the chirp signal 
obtained by setting = 0. 

We use this transmitted chirp signal in (|6]l and apply typical 
demodulation and baseband processing. In particular, the re- 
ceived signal is mixed with the transmitted signal referenced 
to the origin of the scene e"^^'=(*"^'"(^°'+"''(*~^'='(''°)''], 
and then low-pass filtered. If the quadratic phase error is 
ignored lfT4l . we obtain the following signal as the input to 
our algorithm: 



x|l<L 



s(x) e 



-if2fc!(*)[*+*fc-*re/+efci(x,Vx)]vJefci ^-|^2) 



where ftki{t) = ^[<^k - 2ak{t - Tfci(xo))], depends on the 
frequency content of the transmitted waveform. 

The first exponential term depends on the stationary scatters 
only. The second exponential term depends on the moving 
scatters only. For stationary scenes with Vx = 0,Vx, (fT2] i 
represents the 2D Fourier transform of the spatial reflectivity 
function, evaluated at the discrete set of the spatial frequency 
vectors k = [kx,ky]'^ given by 



kfc,i,t = nki{t)eki = ilki(t){ek + ei). 



(13) 



This equation can be used to describe the spatial frequency 
sampling of both the monostatic and the multi-static config- 
uration |9|. Recall that and e; are the unit vectors in the 
direction of the fc-th transmitter and the l-th receiver. For the 
monostatic case these two vectors coincide, i.e. efe = e;. For 
the monostatic case and fixed t, (fT3T l represents an arc of the 
circle of the radius ilki{t) centered the origin of the spatial 
frequency domain, with the length of the arc determined by the 
radial span of the vector e^. Changing t, expands or shrinks 
the circular arcs, leading to the familiar key-hole sampling 
llT4l of the conventional monostatic case. 



For the multi-static case, the situation is a bit more com- 
plicated. For simplicity, assume that flki{t) — ri(t). At fixed 
t and fixed e^, (fT3] l also describes an arc of a circle, this 
time passing through the origin and centered at along 
the direction e^. The length of the arc and its orientation 
are determined by the radial span of the receiver vector e;. 
Changing t results in expansion or shrinkage of the circle 
passing through the origin with its center sliding along the 
direction e^. 

Thus, different covering patters of k-space are possible 
with various sampling strategies in space and time, while 
the resolution is primarily determined by the extent of the 
k-space covering. The range and cross-range resolution of 
both the conventional monostatic SAR system and the multi- 
static, distributed antenna SAR with antenna elements confined 
within the forward cone of Ad < tt/2 are lower bounded by 
the bounding box of the annulus: 



Px > 



2B, 



Pv 



> 



4(,fo + B/2) siii(A6'/2) ' 



(14) 



where Beq = (fo+B /2)-{fa-B /2) cos(Ae'/2). We illustrate 
the k-space sampling patterns used in our experiments in 
Fig. |2] and Fig. [3] Fig. |2|a) and Fig. |2jb) show the k-space 
sampling of the conventional monostatic SAR with bandwidth 
of B = 50MHz, centered at /o = 1.5GHz, over a narrow 
synthetic aperture of A9 = 5 dcg and a wide synthetic aperture 
of A9 — 45 deg, respectively. The k-space covering of the 
narrow-angle monostatic SAR is well approximated with the 
circumscribed rectangle, while this is not the case for the wide- 
angle monostatic SAR. 

A k-space covering similar to the wide-angle mono-static 
case in Fig.|2jb) can be achieved with a multi-static, distributed 
transmit and receive antenna using continuous wave transmis- 
sion. Fig.[3]illustrates the k-space covering for two multi-static 
distributed antenna configurations utilizing continuous wave 
transmission with transmitters and receivers positioned within 
the forward cone of A9 — 45 deg. Fig. Oa) shows a multi- 
static configuration where transmitters sequentially transmit 
a single tone of frequency /o ~ 1.5GHz. Fig. [3jb) shows 
a multi-static configuration where different transmitters send 
out different continuous wave signals within a bandwidth of 
50MHz around the center frequency /q = 1.5GHz. This 
sampling can be achieved either via sequential multi-static 
transmission or via simultaneous MIMO transmission, as such 
signals are easily separated at each receiver Although the 
expected resolutions in the last three cases are similar, it is 
important to notice that the k-space covering of Fig. [3jb) 
for MIMO transmission is achieved in the least amount of 
time, making it the most favorable configuration for imaging 
of moving objects. 

V. Numerical Experiments 

In this section we outline several numerical experiments 
demonstrating the reconstruction capability of the overcom- 
plete dictionary approach for imaging scenes that contain both 
stationary and moving objects. First, we show results for a 
multi-static, distributed antenna sensing configuration for cases 
corresponding to a stationary scene, a scene with moving 
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Fig. 2. Mono-static SAR k-space sampling for B = 50MHz, /o = 1.5GHz with the forward cone centered at Odeg (a) Conventional, narrow-angle SAR 
with Ae = 5deg, (p^ > 2.9m, py > 1.13m) and (b) Wide-angle SAR with AO = 45dcg, (p^ > 0.9m, py > 0.13m) . 
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Fig. 3. Multi-static, distributed antenna SAR k-space samphng with receive and transmit antenna elements positioned within the forward cone of A9 = 45 deg 
centered at Odeg. (a) The case where each transmitter uses continuous waves of frequency /o = 1.5GHz, (px > 1.32m, py > 0.13m) and (b) The case 
where each transmitter sends out continuous waves of different frequencies within the bandwidth oi B = 50MHz centered at /o = 1.5GHz, {px > 0.9m, 
Py > 0.13m) . 



objects whose motion is ignored, a scene with moving objects 
whose motion is expUcitly handled by the overcomplete dictio- 
nary reconstruction algorithm. Next, we show reconstruction 
results for the same cases, but for a wide-angle monostatic 
configuration with the same lower bounds on the range and 
the cross range resolution. Finally, we compare our approach 
with the matched filtering approach described in ID and ifTSll . 

For the multi-static, distributed antenna configuration, all 
transmit and receive elements of the synthetic aperture are 
positioned in the forward cone of AO = 45 deg with its center 
direction aligned with the a;-axis of the coordinate system, as 
illustrated in Fig. [T] We choose a relatively narrow forward 
cone in order to better accommodate the isotropic scatter- 
ing assumption of realistic scatterers. Each antenna element 
transmits a distinct continuous wave signal with a frequency 
randomly chosen within the bandwidth of B = 50MHz around 
the center frequency of /o = 1.5GHz. Similarly, for the 



monostatic case, the angular extent of the synthetic aperture is 
also A6 = 45 deg with the center aspect aligned with the x- 
axis. The transmit waveform is chosen to be the conventional 
chirp signal of _B = 50MHz centered at /o = 1.5GHz. The 
resulting k-space sampling patterns are shown in Fig.|2jb) and 
Fig. [3jb) for the monostatic and the multi-static distributed 
antenna configurations, respectively. 

The scene of interest is of 32 x 32m in size, represented 
by 32 X 128 pixels in the x and y direction respectively, such 
that the spatial cell size is {Ax, Ay) — (1, 0.25)m. The scene 
contains two rigid objects in motion and one rigid stationary 
object, as illustrated in Fig. |4] which displays the ground 
truth images. An object is defined as a set of clustered pixels 
reflecting electromagnetic energy isotropically within the cone 
of 45 deg. There are 20 active scatterers in each object. The 
magnitudes of velocities of the moving objects are 32.5m/s and 
4.7m/s in the direction of tt/6 and tt/6 + tt/2, respectively. 
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Fig. 4. The ground truth of the scene containing one stationary object (upper left corner) and two rigid objects moving at 32.5m/s (lower left comer) and 
4.7m/s (upper right comer): (a) The reflectivity magnitude and (b) The velocity vectors associated with different pixels. 



The reflectivity magnitude of the point scatterers in the scene 
is shown in Fig. |4ja) and their corresponding velocity vectors 
are shown in Fig. lUb). AH measurements are corrupted by 
independent Gaussian receiver noise such that SNR = 20dB, 
with SNR defined as SNR = 20 log J^^^^^^. 

||y-*(V)x||2 

First, we show a set of reconstruction results ifor the multi- 
static, distributed antenna configuration with Nrx = 40 receive 
elements and Nt^ — 10 transmit elements, each sequentially 
transmitting a distinct tone signal with a pulse repetition 
interval of PRI = 2ms. We start by considering what happens 
when there is no motion and when motion is ignored. In 
Fig. |5la) the scene of interest is made completely station- 
ary and our reconstruction is performed with the velocity 
dictionary V containing only the velocity vector vi = 0, 
corresponding to a static scene. In this case, we see that all 
objects are well focused in the reconstruction as would be 
expected. These sparse multi-static stationary scene results 
essentially extend those in presented in |23 | for the case of 
the conventional narrow-angle monostatic SAR. In Fig. |5lb) 
the scene is now made dynamic, as described by Fig. HI but 
scatterers' velocities are ignored in the reconstruction, i.e. the 
dictionary V again contains only the vector vi = 0. We see 
that when the motion of the objects is ignored, the stationary 
object still achieves reasonable focus, while the moving objects 
appear severely blurred. 

Next we demonstrate what happens when we use our over- 
complete dictionary approach to capture the object velocities 
in dynamic scenes. In Fig. [Sjc) we show the reflectivity mag- 
nitude reconstruction when the velocity dictionary contains 
velocities with a magnitude resolution of 3m/s in the range 
[0, 10]m/s and a resolution of Im/s in the range [30,40]m/s. 
The true object velocities are not part of this dictionary. The 
whole dynamic scene is re-focused with scatterer locations 
correctly identified, as illustrated in Fig. EJd). All targets 
appear focused and accurately localized. In Fig. |6] we show 
the corresponding estimated target velocities. All reconstructed 
velocities are correctly estimated within the resolution grid 



error The coarseness of the velocity grid is chosen to avoid 
phase wrapping and further, the phase deviation is minimized 
for shorter CPIs which, in return, are easier to support with 
multi-static and MIMO configurations. 

We now repeat these experiments for the wide-angle mono- 
static configuration. Results are presented in Fig. |7] and Fig.|8] 
showing reconstructed magnitudes and velocities, respectively. 
Recall that each transmitted pulse is now a chirp signal with 
B = 50MHz at fo = 1.5GHz. The number of transmitted 
probes within the angular extent of — 45deg is 40, 
with each pulse return sampled at 10 frequencies. The pulse 
repetition interval is kept as before for consistency (we do 
not worry about the platform velocity required to transverse 
the angular extent within the coherent processing interval). 
The results presented for reconstruction of the dynamic scene 
are thus optimistic. We observe that these reconstructions 
are inferior to those obtained in the multi-static scenario, as 
some features of objects are blurred. The additional diversity 
provided by the multi-static configuration apparently translates 
into improved robustness and quality of the reconstruction. 

Finally, we show reconstruction results for the weighted 
matched filtering approach described in |i9| and Il3l for the 
multi-static, distributed antenna configuration. The reflectivity 
magnitude is evaluated at each pixel for every hypothesized 
velocity, leading to a large space-velocity cube. Here we show 
the maximal value of the reflectivity magnitude of each pixel 
across the velocity direction in this space-velocity cube. 

Fig. |9] shows the results obtained from matched filtering 
when Nrx = 40 receivers and Ntx — 10 transmitters 
are used with ultra-narrow band waveforms, as before. This 
configuration results in 400 observations. Fig. |9^) is the slice 
of the matched filter result at zero velocity, showing what 
happens when velocity is ignored. Fig. |9j)) is the matched 
filter reconstruction of the scene ignoring motion. Fig. |9};) is 
the result obtained by taking the maximum response across the 
space-velocity cube and Fig. |9}l) shows scatter localization, 
by thresholding the result of Fig. |9};) at 0.2 . For the same 
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Fig. 5. Reflectivity magnitude reconstruction with our new, overcomplete dictionary approach for the multi-static, distributed antenna configuration at 
SNR = 20dB with 400 measurements, {Ntx, Nrx) = (10, 40): (a) The reconstruction of the stationary scene assuming no motion, (b) The reconstruction of 
the dynamic scene when velocities are ignored. (V = {vi = 0}). (c) The reconstruction of the dynamic scene with the full overcomplete velocity dictionary, 
(d) The corresponding locations of reflectors in the reconstruction of (c) whose magnitudes are greater than 0.2. 



amount of data, the matched-filter-based reconstructions are 
much worse than those provided by our new overcomplete 
dictionary approach. 

To obtain matched-filter-based reconstructions that have 
similar quality to those produced by the overcomplete dictio- 
nary approach, we can increase the amount of data by using 
Nrx — 40 receivers and Ntx = 10 transmitters, but with the 
transmitters now emitting chirp waveforms of B = 50MHz at 
/o = 1.5GHz with Nf — 30 samples per waveform. Fig. [TO] 
shows the match-filter-based results for this configuration, 
which corresponds to 12,000 observations. Fig. [TOb ) is the 
matched filter reconstruction of the static scene. Fig. [TOb ) is 
the matched filter reconstruction of the scene ignoring motion. 
Fig. [TOb ) is the maximum magnitude reflectively response 
of the matched filter reconstruction and Fig. [TOt l) is the 
thresholded version of this result, showing scatter localization. 
We see that the matched filter reconstruction can come close to 
recovering fine object features, but at the expense of significant 



increase in measurement data relative to the overcomplete 
dictionary approach. 

In Table U we summarize these results by providing the 
per-pixel reflectivity magnitude error for the different cases, 
where this error is defined as E = ||x — x|||/P, with x 
the ground truth reflectivity magnitude, x its estimate and 
P is the number of pixels in the scene. The over-complete 
dictionary results are for 400 observations and the matched- 
filter results are for the 400 observation case and the expanded 
12,000 observation case. The estimates generated by our 
new overcomplete dictionary method are the most accurate, 
producing over a several ten-fold and a 55% reduction in the 
error generated by the matched filter solution which uses the 
same amount of data and 30 times as much data, respectively. 
Further, the spatial diversity of the multi-static, distributed 
antenna configuration provides good reconstructions even in 
the case of narrowband transmission waveforms. 
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Fig. 6. Velocity estimates with our new, overcomplete dictionary approacli for the multi-static, distiibuted antenna configuration at SNR = 20dB. The 
upper right part of the scene with object moving at 4.7m/s: (a) The true velocity field, (b) The corresponding estimate. The lower left part of the scene with 
object moving at 32.5m/s: (c) The true velocity field, (d) The corresponding estimate. 





multi-static, OCD, 
M = 400 


mono-static, OCD, 
M = 400 


multi-static, FBP/MF, 
M = 400 


multi-static, FBP/MF, 
M = 12, 000 


stationary 


0.0015 


0.0018 


0.0452 


0.0040 


dynamic, motion ignored 


0.0104 


0.0129 


0.0646 


0.0093 


dynamic 


0.0036 


0.0090 


0.1378 


0.0085 



TABLE I 

Per pixel reflectivity magnitude error of different reconstruction scenarios with a different number of MEASUREMENTS M. 



VI. Conclusion 

The radar imaging of scenes that contain motion has long 
been an interesting and challenging research topic. We have 
considered multi-static, distributed antenna configurations for 
high-resolution localization of scenes containing both moving 
and stationary scatterers. We have presented an overcomplete 
dictionary inversion approach to simultaneous imaging of sta- 
tionary and moving scatterers. The non-linear, coupled prob- 
lem of joint velocity and reflectivity estimation is effectively 
linearized through introduction of an appropriately defined 
overcomplete velocity dictionary. The resulting optimization 



problem is then approximately solved through convex relax- 
ation. Initial experimental results were presented showing the 
potential of the method for multi-static configurations with 
narrow band transmissions. In contrast to the existing matched- 
filtering approaches that are concerned with the inversion 
of the forward operator only, our overcomplete dictionary 
formulation explicitly and jointly encodes the sparse point 
scattering assumption of both the spatial and the velocity 
dimension, leading to focused imagery with more focused 
object detail. Initial results suggest the new method exhibits 
improved robustness to data loss over existing approaches. 
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Fig. 7. Reflectivity magnitude reconstruction with our approacli for the wide-angle mono-static configuration with SNR = 20dB and 400 measurements: (a) 
The reconstruction of the stationary scene assuming no motion, (b) The reconstruction of the dynamic scene when velocities are ignored, (V = {vi = 0}). 
(c) The reconstruction of the dynamic scene with the full overcomplete velocity dictionary, (d) The con'esponding locations of reflectors in the reconstruction 
of (c) whose magnitudes are greater than 0.2. 
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Fig. 9. Reflectivity magnitude reconstruction of the matclied-filtering/filtered backprojection approach! for tlie multi-static, distributed antenna configuration 
at SNR = 20dB with 400 measurements, {Ntx,Nrx,Nf) = (10,40,1): (a) The reconstruction of the stationary scene assuming no motion, (b) The 
reconstruction of the dynamic scene when velocities are ignored, (c) The maximum reflectivity response in the estimated space-velocity cube for the dynamic 
scene, (d) The conesponding locations of reflectors in the reconstruction of (c) whose magnitudes are greater than 0.2. 
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Fig. 10. Reflectivity magnitude reconstruction of tlie matched-filtering/filtered backprojection approach for the multi-static, distributed antenna configuration 
at SNR = 20dB with 12, 000 measurements, (Ntx, Nrx, Nf) = (10, 40, 30): (a) The reconstruction of the stationary scene assuming no motion, (b) The 
reconstruction of the dynamic scene when velocities are ignored, (c) The maximum reflectivity response in the estimated space-velocity cube for the dynamic 
scene, (d) The conesponding locations of reflectors in the reconstruction of (c) whose magnitudes are greater than 0.2. 



